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Abstract 

We use the 5N formalism to study the trispectrum T( of the primordial curvature perturbation ( when the 
latter is generated by vector field perturbations, considering the tree-level and one-loop contributions. The 
order of magnitude of the level of non-gaussianity in the trispectrum, tnl, is calculated in this scenario and 
related to the order of magnitude of the level of non-gaussianity in the bispectrum, /nl, and the level of 
statistical anisotropy in the power spectrum, . Such consistency relations will put under test this scenario 
against future observations. Comparison with the expected observational bound on tnl from WMAP, for 
generic inflationary models, is done. 

1 Introduction 

Non-gaussianity in the primordial curvature perturbation £ is one of the subjects of more interest in modern 
cosmology, because the non-gaussianity parameters /nl and tnl together with the spectrum amplitude Aq and 
spectral index allow us to discriminate between the different models proposed for the origin of the large-scale 
structure (see for example Refs. [TJ [21 131 H] ) ■ In most of these cosmological models it is assumed that the n-point 
correlators of C are translationally and rotationally invariants. However, since violations of the translational 
(rotational) invariance (i.e. violations of the statistical homogeneity (isotropy)) seem to be present in the data 
[H [SI [71 [HI M Q33 ([H1LI21LI311H11I5]), many researchers have started to build theoretical models that include those 
violations, which could be due to the presence of vector field perturbations [T6l H71 H51 [T9l l20l [211 [22l [23l l24l [25l l26l 
[27J[ii[23[lB[ni[21[Ml[M], spinor field perturbations [351 123 EZ] , or p-form perturbations [HI [Ml HH El EH] , 
due to anisotropic expansion [HI [551 GUI ESH EI3J ESI EH SSI EH HE] or due to an inhomogeneous background 

[M 1301 [Ill- 
Violation of the statistical isotropy is implemented via modifications of the usual definitions of the statistical 

descriptors [16l [49] [50] of the primordial curvature perturbation (. For example, to parametrize the statistical 

anisotropy under the assumption of statistical homogeneity, the power spectrum -P'(k) must include an isotropic 

piece P(k) and an anisotropic piece proportional to the former and exhibiting explicitly the appearance of a 

preferred direction: [5_D] 

p'(k) = p(fc)(i + . 9c (d.k) 2 + ...). (i) 

In the previous expression is a dimensionless parameter, k is the unitary wave-vector, and d is the unitary 
vector along the preferred direction. Some recent papers [Til [HI H31 HH US] claim for the presence of statistical 
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anisotropy in the five-year data from the NASA's WMAP satellite [5T]. In particular, if considering just the 
quadrupolar term of Eq. ^ : 

p c (k) = pf°(fc)(i + . 9c (d.k) 2 ) , (2) 

Ref. [IT] gives g^ ~ 0.290±0.031 which rules out statistical isotropy at more than 9a. Nevertheless, the preferred 
direction lies near the plane of the solar system, which makes the authors of Ref. jTj believe that this effect 
could be due to an unresolved systematic error (among other possible systematic errors which have not been 
demonstrated either to be the source of this statistical anisotropy nor to be completely uncorrelated [H]). Even 
if the result found in Ref. turns out to be due to a systematic error, some forecasted constraints on g^ show 
that the statistical anisotropy subject is worth studying [52]: \g^\ < 0.1 for the NASA's WMAP satellite [51] if 
there is no detection, and \g^\ < 0.02 for the ESA's PLANCK satellite [53 if there is no detection. 

Recent works show that the particular presence of vector fields in the inflationary dynamics may generate 
sizeable levels of non-gaussianity described by /nl [M[[55][56] and tnl [57] ■ As shown in Ref. [56], including vector 
fields allows us to get consistency relations between the order of magnitude of the non-gaussianity parameter 
/nl and the amount of statistical anisotropy in the spectrum g^ . The above studies may transform the violation 
of the statistical isotropy in a decisive tool to discriminate among some of the most usual cosmological models. 
In this letter we use the 5N formalism to calculate the tree-level and one-loop contributions to the trispectrum 
of £ including vector and scalar field perturbations. We then calculate the order of magnitude of the level 
of non-gaussianity in Tf including the one-loop contributions and write down formulas that relate the order of 
magnitude of tnl with the amount of statistical anisotropy in the spectrum, g^, and the order of magnitude of 
the level of non-gaussianity in the bispectrum, /nl ■ Finally, comparison with the expected observational bound 
from WMAP is done. 



2 Trispectrum from vector field perturbations 

The SN formalism [58] [59] [60J [61] [62] extended to include the possible statistical anisotropy in the primordial 
curvature perturbation £ originated from vector field perturbations [3D], provides a powerful tool to calculate £ 
and its statistical descriptors. Assu min g an inflationary dynamics dominated by just one scalar field 4> and one 
vector field A, £ is expressed as |30j#*l 



C(x) = flV(0(x), Mx), t) = N^S</> + N\5Ai + ^N^(6<f>) 2 + N^SA, + ^N^ A 6A t 6A 3 , (3) 

where 

dN . ON Ar d 2 N ~ d 2 N . d 2 N 

and i denotes the spatial indices running from 1 to 3. Now, we define the power spectrum and the trispectrum 
for the primordial curvature perturbation, through the Fourier modes of £, as: 

2tt 2 

(C(ki)C(k 2 )) = (27r) 3 <5(k 1+ k 2 )P c (k) = (2 7 r) 3 5(k 1 + k 2 )— P c (k) , (5) 
(C(k,)C(k 2 )C(k,)C(k 4 )) = (2 7 r) 3 ^(k 1 +k 2 +k3 + k 4 )r c (k 1 ,k 2 ,k3,k 4 ) 

= (27r) 3 5(k 1 +k 2 +k 3 +k 4 )^pJ^- 1 -^rc(k 1 ,k 2 ,k3,k4). (6) 

As shown in Ref. [3D] , the tree-level contribution to the spectrum has the form of Eq. @ . This is simply obtained 
by using Eqs. ([3]) and ([5]). Assuming again only tree- level contributions, the level of non-gaussianity /nl in the 
bispectrum was calculated in Ref. [54]. The same calculation was performed in Ref. [56] but this time 
including also one-loop contributions and considering them to be dominant over the tree-level terms. In both 
works, and /nl were shown to exhibit anisotropic contributions coming from the vector field perturbation. 
In this letter we show that it is possible to obtain an analogous expression for the level of non-gaussianity ttstl 



# lr This expression corrects Eq. (3.14) of Ref. 1301 . and Eq. (3) of Ref. |54| . where a factor 2 in the fourth term of the expansion 



is missme 
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in the trispectrum Tq. To do it, we first need to calculate t he ex pressions for Pq and Tq, defined in Eqs. ([5]) and 



©• Considering contributions up to one-loop order, we find# 2 



7f ee (k) = NlVsm + N^NiTiiW 

= Np>s4>(k) + N 2 A V+(k) + (N A ■ k) 2 V+(k) (r long - 1) 



(7) 



V, 



1 — loop 



(k) 



d 3 p k 3 



47r|k + p| 3 p 3 



-N^Vs^lk + p\)Vst(p) + Ni A Ni A V s<t> (\k + p|)7^(p) 



-2^M^(k + p)^(p) 



(8) 



7? roc (k 1 ,k 2 ,k 3 ,k 4 ) 



= N^N^[Vs^k 2 )V8^{ki)Vs^{\ki + k 2 |) + 11 perm.] 

+ N A N A N AA N AA [7- fe (k 2 )7^ m (k 4 )7i„(k 1 + k 2 ) + 11 perm 

+ N%Njy,N s M [psiWPBtikJTijQiLi + k 2 ) + 11 perm 

+ N A N A N^N A ^Ti k (k 2 )Tjifa)P S 4\ki + k 2 |) + 11 perm 

+ N^N^NXN^ ^(fe)75j(k4)P^(|ki + k 2 |) + 23 perm 

+ Nj) N A N A( p N AA IPs^TiMTjiQn + k 2 ) + 23 perm 



(9) 



where 
and 



T{ A loop (ki,k 2 ,k 3 ,k4 



Ajij njkl N mn N op 



d 3 p fcifc||k 3 + k 4 



47rp 3 |ki - p| 3 |k 3 + p| 3 |k 3 + k 4 + p| 3 
x7I m (p)7}fe(ki - p)7; ip (k 3 + p)7Io(k 3 + k 4 + p) , 



Ty(k) = T^(k)P + (k) + iT9f d (k)P.(k) + T^(k)P long (k) . 



r?y™(k) = % - hk 



3 ' 



T ij (k) = eyfcfcfc , 



T>° n S(k) = fc . fc 



(10) 

(11) 
(12) 



Eq. was written in the form of Eq. (0) with d = N^, being a vector with magnitude N A = -J ' N A N l Al 
and ri on g = Piong/P+, where Pi ong is the power spectrum of the longitudinal component, and P+ and "P- are 
the parity conserving and violating power spectra defined by 



P ± = -(P R ±P L ) , 



(13) 



with Pr and Pl denoting the power spectra of the transverse components with right-handed and left-handed 
polarisations [30]. Eq. (fTU)) only includes terms coming from vector field perturbations; this is because the 
complete expression (including the scalar and the mixed terms) is too large and in the current paper we are 
assuming that the contributions to T> coming only from vector fields dominate over all the other contributions. 



3 Vector field contributions to the statistical descriptors 

When statistical anisotropy is assumed, there is an important restriction from observation: one related to the 
amount of statistical anisotropy present in the spectrum, which is given by the parameter g^ in Eq. ([2]). Recent 
studies of the data coming from the WMAP experiment, set an upper bound over g^ < 0.383 [IT] . The latter 
observational constraint is fully satisfied when we assume that the contributions coming from vector fields in 

# 2 Eq. {5J corrects a mistake in Eq. (4.12) of Ref. 1301 where the infinitesimal volume element cfip was incorrectly expressed in 
terms of dp. 
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Eqs. (0 and ([5]) are smaller than those coming from scalar fields. That means that the first term in Eq. ([7]) 
dominates over all the other terms, even those coming from one-loop contributions. 

In our study we will assume that the terms coming only from the vector field dominate over those coming 
from the mixed terms and from the scalar fields only, except for the case of the tree-level spectrum, where we 



will assume that the scalar term is the dominant on qj_f| . Of course, for an actual realisation of this scenario, we 



need to show that such constraints are fully satisfied. From the above assumptions it follows that: 

7> c trcc (k) = p c *"(k)+p c f°Q4, (w) 

pl-loap (k ) = 7 , c ^l°oP( k ) j ( 15 ) 

^(ki.ka.ka.k*) = r c t ; oc (k 1 ,k 2 ,k3,k 4 ), (16) 
^"^(kx.ka.ka.k*) = TcV loop (ki,k 2 ,k 3 ,k 4 ), (17) 

where the subscripts ^ and C, A mean scalar field or vector field contributions to £. The above expressions lead 
us to two different possibilities that let us study and probably get a high level of non- gaussianity: 

• Vector field spectrum (V^ A ) and trispectrum (7^ A ) dominated by the tree-level terms. 

• Vector field spectrum (V^ A ) and trispectrum {Tq a ) dominated by the one-loop contributions. 

Other possibilities are not viable because it is impossible to satisfy simultaneously that the vector field spectrum 
(V^ A ) is dominated by the tree-level terms and the trispectrum (Tq a ) is dominated by the one-loop contributions, 
or the vector field spectrum (Vc . ) is dominated by the one-loop contributions and the trispectrum (Tq a ) is 
dominated by the tree-level terma# 4 l. This is perhaps related to the fact that we have taken into account only 
one vector field. Such a conclusion may be relaxed if we take into account more than one vector field, as 
analogously happens in the scalar multi- field case [531 El] ■ 

In order to study the above possibilities, we need to estimate the integrals coming from loop contributions. 
From Eqs. CDS, dHJ), and (JTrJ) the integrals to solve are: 

n^ lo ° P (k) = J«/^a^ + P)^(p), (18) 

,-pi-ioop/, , , , n _ AT ij AT u i\rmn AT o P f d 3 p fcffcf |k 3 + k 4 | 3 

T CA frMMM) - N AA N AA N AA N AA J _ ^ + ^ - - - ^ x 

x7I m (p)7}/t(k 1 - p)7; ip (k 3 + p)7Io(k 3 + k 4 + p) . (19) 

The above integrals cannot be done analytically, but they can be estimated using the same technique shown in 
Appendix El it is found that the integrals are proportional to hx(kL) (where L is the box size) if the spectrum 
is scale invariant. Following it, we find from Eqs. ([T5]) and flTT)l) : 

^V OOP (k) = \N i l A Nf A {2V + + V, ms )8 ik r j mHkL), (20) 

Tl A loop Qti, k 2 , k 3 , k 4 ) = Nl A Nf A N%2N° A p A \n(kL)(2T+ + V long )5 tm [7^ fc (k 1 )r„ P (k 3 )7I (k 4 + k 3 )] . (21) 

Observations are available within the observable universe and, except for the low multipoles of the CMB, all 
observations probe scales k 3> Hq. To handle them, one should choose the box size as L — H^ 1 [65]. A smaller 
choice would throw away some of the data while a bigger choice would make the spatial averages unobservable. 
Low multipoles £ of the CMB anisotropy explore scales of order H^ 1 /£ not very much smaller than Hq 1 . To 
handle them one has to take L bigger than Hq . For most purposes, one should use a box, such that \n(LHo) 
is just a few (i.e. not exponentially large) [66j ETJ [68] . When comparing the loop contribution with observation 
one should normally set L = Hq 1 , except for the low CMB multipoles where one should choose L 3> -f^o" 1 with 
ln(fci) ~ 1. With the choice L — Hq , ln(fcL) ~ 5 for the scales explored by the CMB multipoles with I ~ 100, 
while ln(fci) ~ 10 for the scales explored by galaxy surveys. Since we are interested in giving orders of magnitude 
and simple mathematical expressions, we will set ln(fcL) ~ 1 without loss of generality. 



#3 Thc 

power spectrum must be dominated by the tree-level terms. Otherwise there would be too much scale dependence in 
conflict with the current observational limit on n^. 
# 4 See the relevant discussion regarding the vector field bispectrum B^ A in Ref. |56| . 
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4 Calculation of the non-gaussianity parameter tnl 



The non-gaussianity parameter tnl is defined by [69) : 

= 2T c (k 1 ,k 2 ,k 3 ,k 4 ) 

TNL [n(ki)n( k 2)n( k i + k 4) + 23 perm.] ' 1 ' 

Remember that the isotropic contribution in the Eq. ([2]) is always dominant compared to the anisotropic one so 
that we may write in the above expression only the isotropic part of the spectrum V^ lso (k): 

= 27^(ki,k 2 ,k 3 ,k 4 ) 

TNL [P c iso (fc 1 )P < iso (fc 2 )P c iso (|k 1 +k 4 |) + 23perm.] ' 

Using the above expression, we will estimate the possible amount of non-gaussianity generated by the anisotropic 
part of the primordial curvature perturbation, taking into account different possibilities and assuming that the 
non-gaussianity is produced solely by vector field perturbations. 

4.1 Vector field spectrum (V^ A ) and trispectrum (7^ 4 ) dominated by the tree-level 
terms 

In this first case, we assume that the trispectrum is dominated by vector field perturbations and that the higher 
order terms in the 5N expansion in Eq. ([3]) involving the vector field are sub-dominant against the first-order 
term: N A SAi ^> N AA SAi6Aj. The latter implies that both the spectrum and the trispectrum are dominated by 
the tree-level terms, i.e. V c ^ cc > V c 1 A loop and T Ca cc > T c Jf loop . Thus, we have from Eq. (23]): 

= 2rc t j cc (ki,k 2 ,k3,k 4 ) 

NL [V c iso (k 1 )T c iso (k 2 )T c iso (\k 1 + k A \) + 23perm.] , 

which, in view of Eqs. © and (|T6"|) . looks like: 

7^(k 2 )75 m (k 4 )7I„(ki + k 2 ) + 11 perm. 



2N\N J A Nf A N^2 



[n iS °(^i)n iS0 (fc2)n is °(|ki + k 4 |) + 23 perm.] 



(25) 



We will just consider here the order of magnitude of tnl- Therefore, we will ignore the specific k dependence 
of Tij- Instead, we will assume that V\ ng, V+, and V- are all of the same order of magnitude, which is a good 
approximation for some specific actions (see for instance Ref. [30]), and take advantage of the fact that the 
spectrum is almost scale invariant |70) . Thus, after getting rid of all the k dependences, the order of magnitude 
of tnl looks like: 

T NL ~ A A AA / 2 6) 

where V a = 2V+ + Piong- Employing our assumption that NaSA > NaaSA 2 , and since the root mean squared 
value for the vector field perturbation 5 A is \jPa~, the contribution of the vector field to £ is given by Qa ~ 
\/Vc A ~ Na\ZVa- An upper bound for tnl is therefore given by: 

tnl < VCA ■ (27) 

CP c ISO )3 

Since the order of magnitude of g^ is 'P^/'P^ 1 ' 30 , under the assumptions made above we get: 

where ("P^ ) 1 / 2 ~ 5 x 1CP 5 [70] has been used. Eq. (|28p gives an upper bound for the level of non-gaussianity 
tnl i n terms of the level of statistical anisotropy in the power spectrum g^ when the former is generated by the 
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anisotropic contribution to the curvature perturbation. Com pari ng with the expected observational limit on tnl 



coming from future WMAP data releases, tnl ~ 2 x 10 4 [7lJ# 5 l. we conclude that in this scenario a large level 
of non-gaussianity in the trispectrum of £ is possible, leaving some room for ruling out this scenario if the 
current expected observational limit is overtaken. 

As an example of this scenario, we apply the previous results to a specific model, e.g. the vector curvaton 
model [171QI1II2], where the ^-derivatives are [54]: 

Na = ^r, (29) 

Naa = jV r > ( 3 °) 

where A = |A| is the value of vector field just before the vector curvaton field decays and the parameter r is the 
ratio between the energy density of the vector curvaton field and the total energy density of the Universe just 
before the vector curvaton decay. 

First, we check if the conditions under which the vector field spectrum and trispectrum are always dominated 
by the tree- level terms are fully satisfied. From Eqs. ([7]), ©, ((2"0|) and (f2"Tj) our constraint leads to: 

V A N\ » V 2 A N\ A , (31) 

KKKa » Kn\ a , (32) 

which mean that the if the vector field spectrum is dominated by the tree-level terms so is the vector field 
trispectrum. An analogous situation happens when the vector field spectrum is dominated by the one-loop 
terms: the vector field trispectrum is also dominated by this kind of terms. As a result, it is impossible that 
simultaneously the vector field spectrum is dominated by the tree-level (one-loop) terms and the vector field 
trispectrum is dominated by the one-loop (tree-level) terms. Following Eq. (|3ip . we get: 



which, in view of (a ~ y/^CA ~ NaVPa and Eqs. ([29]) and (|30]) . reduces to: 

r > 2.25 x l(r 4 ^ /2 . (34) 

This lower bound on the r parameter has to be considered when building a realistic particle physics model of 
the vector curvaton scenario. 

Second, looking at Eq. (l26l) . we obtain the level of non-gaussianity tnl for this scenario: 



2 x 10" 2 

r NL — n 



m 3 ■ < 35 » 



This is a consistency relation between tnl , .<?£ j an d r which will help when confronting the specific vector curvaton 
realisation against observation. Indeed, a similar consistency relation between /nl and was derived for this 
scenario in Ref. [55] : 

4.5 x IP' 2 ( 9i \2 
/nl^ (^J • (36) 

Thus, in the framework of the vector curvaton scenario, the levels of non-gaussianity /nl and tnl are related to 
each other via the r parameter in this way: 

2-1 ,3/2 /on ,-, 

T NL ^ -J72/NL ) K 67 ) 

in contrast to the standard result 

r NL = t^t/nl > ( 38 ) 

for the scalar field case (including the scalar curvaton scenario) found in Ref. [73]. 



# 5 The trispectrum in this scenario might be either of the local, equilateral, or orthogonal type. We are not interested in this letter 
on the shape of the non-gaussianity but on its order of magnitude. Being that the case, comparing with the expected bound on the 
local tnl |71| makes no sensible difference under the assumption that the expected bounds on the equilateral and orthogonal tjml 
are of the same order of magnitude, as analogously happens in the /nl case for single-field inflation 72 . 
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4.2 Vector field spectrum (V^ A ) and trispectrum (7f x ) dominated by the one-loop 
contributions 

From Eqs. and (J33J) we get 

„ N l j A NX l A N^N / A HkL)(2V + + P 1 on g )^ m [^ fc (k 1 )r rap (k3)7I (|k 4 + k 3 |)] 
TNL [P c iso (A ;i )P c iso (fc 2 )P c iso (|k 1 +k4|) + 23perm.] 

Assuming again that 7-Wg, V+, and T 3 - are all of the same order of magnitude, and that the spectrum is 
scale invariant, we end up with: 

-p 4 /V 4 

TNL ~ U^AA (40) 

Performing a similar analysis as done in the previous subsection, but this time taking into account that the 
vector field spectrum is dominated by the one- loop contribution and therefore £4 ~ "v^Ca ~ ^AaPa-, we arrive 
at: 



2 



^ 8 x 10 6 



(V c iSO ) 3 

The above result gives a relation between the non-gaussianity parameter tnl arid the level of statistical anisotropy 
in the power spectrum g^. 

Now, we call a similar result that we found for the non-gaussianity parameter /nl in Ref. |56] . that is: 

/n L ~10 3 (^) 3/2 . (42) 

By combining Eqs. fill) and p2|) we get: 

t nl ~ 8 x 10 2 /nl ; (43) 
which gives a consistency relation between the non-gaussianity parameters /nl and tnl for this particular 
scenario. The consistency relations in Eqs. (|41[) . (|42|). and (|43| will put under test this scenario against future 
observations. In particular, the consistency relation in Eq. ("43)) differs significantly from those obtained when £ 
is generated only by scalar fields (see e.g. Eq. (|38| and Ref. [73]). 

Again when we apply our result to the vector curvaton scenario, we get from Eqs. ([7]), ©, (f2"0)) . (|2T|) . ([29]) 
and (J3UJ) : 

r < 2.25 x 1(T 4 ^ /2 , (44) 

which is an upper bound on the r parameter that must be considered when building a realistic particle physics 
model of the vector curvaton scenario. 



5 Conclusions 

We have studied in this letter the order of magnitude of the level of non-gaussianity tnl in the trispectrum Tf 
when statistical anisotropy is generated by the presence of one vector field. We have shown that it is possible to 
get an upper bound on the order of magnitude of tnl if we assume that the tree-level contributions dominate 
over all higher order terms in both the vector field spectrum (V^ A ) and the trispectrum (7^ A ); this bound is 
given in Eq. (|2"8")l . We also show that it is possible to get a high level of non-gaussianity tnl, easily exceeding 
the expected observational bound from WMAP, if we assume that the one-loop contributions dominate over the 
tree- level terms in both the vector field spectrum (V( A ) and the trispectrum (T( A ) . tnl is given in this case by 
Eq. (l4Tj) . where we may see that there is a consistency relation between the order of magnitude of tnl and the 
amount of statistical anisotropy in the spectrum g^. Two other consistency relations are given by Eqs. (|42|) and 
(j43| , this time relating the order of magnitude of the non-gaussianity parameter /nl in the bispectrum Bq with 
the amount of statistical anisotropy and the order of magnitude of the level of non-gaussianity tnl in the 
trispectrum Tq. Such consistency relations let us fix two of the three parameters by knowing about the other 
one, i.e. if the non-gaussianity in the bispectrum (or trispectrum) is detected and our scenario is appropriate, 
the amount of statistical anisotropy in the power spectrum and the order of magnitude of the non-gaussianity 
parameter t N l (or /nl) must have specific values, which are given by Eqs. (|4"2")) (or ([H])) and ([4*3")) . A similar 
conclusion is reached if the statistical anisotropy in the power spectrum is detected before the non-gaussianity 
in the bispectrum or the trispectrum is. 
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A One-loop integral for Pq 

We sketch in this appendix the mathematical procedure to estimate the integrals in Eqs. © and ([TO)) . We only 
work one integral since the other ones are estimated in a similar way. 
The one-loop contribution to the spectrum is: 



d 3 p k 3 

n 

~iV^< A ^(k + p)^(p) 



(45) 



As we can see, the total contribution to "pi~ loo P corresponds to three integrals, each one having two singularities: 
one in p = and the other one in p = — k. If the fields spectra are scale invariant, the first integral may be 
written as: 

"^w-s^/i^&F- (46) 

so the actual integral to estimate is: 

f d 3 p k 3 

1 = J L -, |k + p|3p3 • ( 4? ) 

This integral is logarithmically divergent at the zeros in the denominator, but there is a cutoff at k = L^ 1 . The 
subscript L" 1 indicates that the integrand is set equal to zero in a sphere of radius L~ x around each singularity, 
and the discussion makes sense only for L _1 <C k <C k max . If we consider the infrared divergences, that means 
p -C k, we may write: 

/ = / ~47rln(feL). (48) 

7x,-i p 

To calculate the contribution coming from the other singularity we can make the substitution q = k + p. After 
evaluating this latter integral, we find that the contribution is again 47rln(fcL). The integral in Eq. (|47|) may be 
finally estimated by adding the contributions of the two singularities, that means: 

/ = / ,/ P l s =87rln(fcL). (49) 

J |k + phV 

More details to evaluate these integrals may be found in Refs. [STJ IM1 174] . 

The technique to evaluate this kind of integrals when considering vector fields is the same, although the 
procedure is algebraically more tedious. Nevertheless, one can finally arrive to the same conclusion. A more 
detailed discussion about this issue will be found in a forthcoming publication [75] . 
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